Abstract 



We investigate the dynamics of spin-nonequilibrium electron systems 
for the case when normal electron collisions prevail over the other scat- 
tering processes and the "hydrodynamic flow" regime is realized. The 
hydrodynamic equations for the electron liquid have been obtained and 
analyzed. We demonstrate that oscillations of the spin polarization are 
possible in a conducting ring with inhomogeneous magnetic properties. 
These low-decay oscillations are accompanied by the oscillations of the 
drift current in the ring. We demonstrate also that the spin polarization 
of the electron density may be revealed via the voltage between the ends 
of the open circuit with an inhomogeneous spin polarization. The effect 
may be observed both in the hydrodynamic and diffusive regimes. 
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Introduction. Last years, generation of nonequilibrium spin-polarized elec- 
tron states in non-magnetic conductors and manipulation with them were in the 
focus of attention (see, e.g. [1]). A number of promising methods and schemes 
was proposed and discussed; the most of them operate under the conditions 
when electron-electron interactions are neglected. Effects of electron-electron 
scattering in the spin dynamics were considered in Refs. 2-4, but it was not for 
the case of hydrodynamic electron flow ^ , when the electron-electron scatter- 
ing dominates and the electron system may be considered as a liquid with its' 
inherent effects. 

Meanwhile, a hydrodynamic flow regime is quite real in a low-dimensional 
electron gas in heterostructures [6 (as well as, in electron systems over the liquid 
helium surface [7J). The main condition is the following. The transport electron 
mean free path should be determined by the "normal" collisions that conserve 
the total momentum of the system of interacting particles. It may be electron- 
electron collisions (when the Umklapp processes are absent due to small sizes of 
the Fermi surfaces of low dimensional conductors) or electron-phonon collisions 
(when phonons are tightly coupled to the electron system). In other words, the 
condition 1^ << ly should be satisfied, where In is the mean free path (m.f.p.) 
for the normal collisions, ly is m.f.p. for bulk collisions that do not conserve 
the quasi-momentum. 

A hydrodynamic electron flow was observed experimentally in a high- 
mobility electron gas in (Al,Ga)As heterostructures (SI (see, also Ref .[9 ) at 
1,5 - 20 K. For that range of the electron temperatures the electron-electron 
m.f.p., lee, is much less than the electron m.f.p. determined by the collisions 
with imperfections of the heterostructure, h. Beside that, in the experiment [5] 
the phonon temperature was much lower than the electron temperature and col- 
lisions with phonons were unessential. The hydrodynamic regime breaks under 
increase of the temperature of the sample, because the electron-phonon m.f.p., 
lep, becomes less than lee, and phonons can remove effectively a momentum 
from the low-dimensional electron system being in a good acoustic contact with 
the surrounding media. If, in the experiment, it is possible to provide a weak 
contact with the surrounding media (that results in the conservation of the total 
momentum of the electron-phonon system) , then the increasing of the tempera- 
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ture favors the hydrodynamic regime due to the decreasing of the normal m.f.p., 

In the hydrodynamic flow regime, the state of the electron liquid is charac- 
terized by the velocity u(r) and its density p{r). Naturally, when our electron 
system is spin-polarized, the densities of the spin components differ from each 
other. In this case the electron liquid is a two-component mixture and the 
density variables have the spin indexes, i.e. Pcr(r). Meanwhile, in the leading 
approximation, the velocity u(r) is the same for the all spin components. The 
reason is that frequent collisions between electrons with different spins form a 
common drift of the electron system. Moreover, we have to regard the electron 
liquid as incompressible when the geometric size of a conductor is larger than 
the electron screening radius (which is comparatively small in metals and het- 
erostructures) and the characteristic frequencies of the considered processes are 
less than the plasma frequency. We would like to emphasize that the incom- 
pressibility means here that the nonequilibrium addition to the electron density 
is vanishing (nevertheless, the equilibrium density may be inhomogeneous due 
to the reasons that are discussed below). Obviously, the incompressibility and 
the equal velocities of the spin components lead to the following simple fact. 
The total current through the channel cross-section, /, does not depend on the 
coordinate along the channel and it is distributed between the spin components 
in proportion to their densities. Thus, the current characteristics of the system 
are determined by the value of /. 

As we demonstrate below, in the conducting ring with inhomogeneous mag- 
netic properties the value of the total current undergoes low-decay harmonic 
oscillations which frequency depends on the characteristics of the inhomogene- 
ity. The nature of these oscillations is the following. The drift in the spin- 
inhomogeneous ring causes appearing of a non-equilibrium spin polarization, i.e. 
an accumulation of the non-equilibrium densities of the spin-up and spin-down 
components occurs (while the total density is conserved). The accumulation 
exists until the moment when the drift is stopped due to the interaction of non- 
equilibrium electrons with a field that induces the inhomogeneity of the electron 
spectrum. But, electrons have the inertial masses and the process will evolve 
back. We named this oscillation process the "spin pendulum". 

Note, the existence of the well-known hydrodynamic waves, which frequen- 
cies are less than the plasma frequency, is impossible here due to the Coulomb 
interaction. Consequently, " spin pendulum-like" oscillations are the only possi- 
ble oscillations of the system in this case. 

The interesting spin-electrical effect related to the non-equilibrium spin po- 
larization may be observed not only in the ring but in the open circuit as well. 
The skewed spin polarization causes a voltage between the ends of the open 
circuit and it allows easily registering of the spin polarization. Note, the spin- 
electrical effect in magnetoelectric materials in the equilibrium inhomogeneous 
state was discussed early in [TD]. 

Note that the all aforementioned effects are possible both in the magnetic 
conductors, where densities of the spin components differs from each other ini- 
tially (e.g., due to the exchange interaction), and in non-magnetic materials, 
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where either spin separation is due to the Rashba effect [TT], or it is caused 
by an external magnetic field. It is important that non-magnetic materials get 
magnetic properties due to the appearing of the induced non-equilibrium spin 
distribution. Consequently, the aforementioned steady-state and dynamic ef- 
fects exist in them but it is the second-order effects in the spin non-equilibrium. 
In the case of a two-dimensional electron gas in heterostructures, the variation 
of magnetic properties can be induced by a non-uniform gate voltage applied, 
which affects the Rashba effect [12] ; by a variation of the external magnetic field 
or by a space-dependent spin injection. 

Spin hydrodynamics A time-dependent distribution function /(r, p, t) for 
electrons at the position r and with the momentum p obeys the Boltzmann 
transport equation 

dt dr dr dp ^' 

The hydrodynamic equations can be derived from the Boltzmann transport 
equation by the method used in Refs. [131 El) taking into account that the 
function /(r,p,t) depends also on the spin index a, which corresponds to the 
different spin components. It is assumed in Eq.([T]) that the energy spectrum 
of the current carriers depends on the coordinates, momentum and spin index: 
e = e(r,p, cr). The electrical potential ip appears due to the space-dependent 
non-equilibrium spin density. J(/) is the scattering term and v = de/dp is the 
electron velocity. 

Let us assume that momentum dissipation is vanishing and the normal scat- 
tering processes dominate, and the following conditions are fulfilled: In/L « 1 
and uitn << 1. Here L is the characteristic geometric size of the system, ui is the 
characteristic frequency of the oscillation processes in the system, tn = /vp is 
the relaxation time corresponding to the normal collisions, and vp is the Fermi 
velocity. In this case we may expand Eq.(l) in series in these small parameters. 
In the leading approximation we obtain the following equation 

JnU'^'^) = 0. (2) 

Here Jjv is the part of the collision term which corresponds to the normal 
collisions. The solution of Eq.([2]) is the quasi-equilibrium drift distribution 

/(O) = n{e, - 5fi, - up), n(z) = {e^^-^^)/^ + (3) 

where the drift velocity u and nonequilibrium additions to the chemical poten- 
tial Sfia depend on the coordinates r and time t (we neglect here the temperature 
variations). Following the method used in Rcfs. [131 [Hj we analyzed the con- 
ditions of solvability up to the next two orders of series expansions and have 
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obtained the following equations for Sfi^j^u, and ip 

86 Pa 



dt 



du 



dt ^ 



Per = PaiPa)- 

Here po- = pqct + "^Pt is the spin-dependent current density, and tq^ is the 
equilibrium spin density. (Note that 5pa is the function on 5^,^, and, in the linear 
approximation, 5pa — lia^p-m where lia is the spin-dependent density of states 
on the Fermi surface.) mp = (x)^ "12 I ^^P^ P — Pa where r is 

the dimensionality of the electron system. D is the diffusion operator; F is the 
spin-flip operator; the operator V relates to the viscosity and the operator I] 
corresponds to the scattering with momentum dissipation (see, e.g. |131 114|). 
Note that the specific forms of these operators are not important for us: the 
corresponding terms in Eqs.(|][]) are small enough and we use these terms only 
for the estimates. 

Meanwhile, Eqs. ([]■[]) have a clear physical meaning. Equation ([]) describes 
the law of conservation of the number of colliding particles, while equation ([| 
corresponds to the conservation of the total momentum. Equation ([]) is the 
result of incompressibility of the electron liquid with the Coulomb interaction. 
(To derive equations (4-5) in the form presented, we assume that electrons have 
a spherical Fermi surface.) Eqs. (|][]) are written in the linear approximation in 
the drift velocity u, but the terms that contain bp„^ 5pa are given in the explicit 
form under the assumption u — > 0. This approximation is sufficient here since 
the non-linear effects will be considered for the case when the common drift is 
absent. 

Neglecting both the viscosity of the electron liquid and momentum dissipa- 
tion we rewrite Eq.([]) in the following form {5n — n{e — Sfi^) ~ ^(s)) 

mp— + VP + epVip + }2 —^'^P = 0' (^) 



pvfe dP 



Note that, integrating Eq.Q by parts in the momentum space, we obtain Eq.([]). 
Equation Q is a natural generalization of the Euler equation T5] for the mixture 
of liquids interacting with an external electrical field and with the field that in- 
duces the inhomogeneity of the electron spectrum. P{5pa) is a non-equilibrium 
addition to the electron pressure. 

In the case of a conducting ring, the electron pressure does not cause the net 
effect on the electrons, this is because § VPdx = 0(if our conductor is homoge- 
neous, pressure can be compensated by an electrical field, V</?). Meanwhile, the 
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Figure 1: Equilibrium distribution of the spin density po-f along the ring of length 
L (Fig. la). When electrons drift into the right (see Fig. lb) the non-equilibrium 
spin-up potential S^^ appears in the regions where the spin-up density varies. 
The total force is a sum of all forces acting on electrons and it is directed to the 
left. Consequently, it leads to the change of the drift direction. 



force VEct may cause a common drift in the ring as it acts on the electrons which 
are non-equilibrium on the spin. Figure 1 illustrates the physical picture arising 
in the conducting ring with the inhomogeneity of the equilibrium spin density 
and explains the nature of the total force ating on the spin-nonequilibrium dis- 
tribution caused by the electron drift along the inhomogeneous ring. 

Spin pendulum. Let us discuss the oscillation process described by Eqs. 
m which is related to the current / = /qc*"*, i.e. a spin pendulum. In the linear 
approximation we may solve the problem completely and find the frequency of 
the oscillations 
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(7) 



Here x is the coordinate along the conducting channel, " up arrow" corresponds 
to the one of spin components; s is the cross-section of the channel which is 
assumed a;-independent; strike denotes differentiation with respect to x, L is 
the length of the conducting channel (i.e. the ring length, in our case). 

As is evident from Eqs.(6][7l, the oscillations exist only when the rela- 
tive concentration of the spin components varies along the channel. The fre- 



6 



quency of oscillations is proportional to the level of the magnetic inhomogeneity, 

a = {P-y-niax - P<Tmj«)/P> viz. 

Lo Ki avF{LLi) ""^j (8) 

where Li is the characteristic scale of the spectrum inhomogeneity. 

Initially, the oscillations of the spin-pendulum could be exited by the " mag- 
netic push" when an external magnetic field is switched on (much faster than 
Lo~^) and an initial current / = es^/cra is induced (where $ is the magnetic 
flux through the ring). Naturally, there are dissipation processes in the hydro- 
dynamic flow, which we did not take into account. They lead to the damping 
of the spin-pendulum oscillations. It is easy to estimate the decay time 



/ e^psR D I \ ,„. 



where R is the electrical resistance of the homogeneous ring. The second term 
in the r.h.s. of Eq.(|9| is related to the diffusion of the nonequilibrium spins 
through the inhomogeneous region. In the hydrodynamic regime the diffusion 
coefficient is determined by the normal collisions, i.e. D « I^vf- The third 
term in the r.h.s. of Eq.Q corresponds the spin relaxation which is due to the 
spin-flip with the characteristic time Tgf. Note, as it follows from Eq.([9|, the 
"space-sharp" inhomogeneity (e.g., a sharp interface boundary between different 
magnetic materials) causes a very fast decay of the oscillations. 

It seems interesting that there exists a possibility to keep the amplitude 
of the oscillations constant (like in an ordinary pendulum clock) owing to the 
magnetic connection of the ring with an energy-pumping cell. 

Spin-electrical effect. Let us discuss now a spin-electrical effect in an 
open-ended homogeneous conductor. At u = we obtain from Eqs. (|4f - 
the following potential difference between the open ends 

^(L)-^(0) = -[P(0)-P(i)], (10) 
P 

P{-) - f^t - ^Vmt + fn, - ^g) + ... (11) 



Eq. (Ill yields the electron pressure as the expansion in series on the (5/i. 
The second order term is written in the form for the case of non-magnetic 
materials. All the quantities in the coefficients of Sjj,^ and [5^^)"^ are understood 
as the equilibrium state quantities. It is obvious from Eq.(|lO|) that the difference 
between the non-equilibrium spin pressures at the ends of the circuit induces 
an electrical voltage that can be measured experimentally. The voltage exists 
during the lifetime which is determined either by the spin relaxation due to 
the spin-flip processes or by the diffusion equalization of the spin concentration 
along the circuit, « {t~^ + D / L'^)^^ . The electrical charge q « Aipr^ flows 
through the meter during this time. Let us discuss the case when the spin 
polarization is induced in a non-magnetic conductor due to the current spin 
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injection from the magnetic material. Then, a voltmeter connected to the ends 
of the non-magnetic conductor will fix the electrical charge q and indicate the 
non-equilibrium spin density be existed in the non-magnetic conductor. 

We should note here that this "spin-electrical" effect is not a specifically hy- 
drodynamic effect. In the diffusion regime, when collisions that do not conserve 
momenta are dominated over normal collisions, we may write the continuity 
equations in the following form 

diva,(V((5Ai, + e(p)) =0, (12) 

where aa- is the contribution to the conductivity of the corresponding spin com- 
ponent, a = ^cFa-- Consequently, taking into account Eq.([l) we obtain the 

following equation for the case of homogeneous conductor 



511+ 9(Tt\ 



,2 



(13) 



Here, the second term in the r.h.s. of Eq.(19) is written for the case of non- 
magnetic conductors only. 

Conclusion. In summary, the low-decay oscillations of the spin polariza- 
tion with the frequency w accompanied by the oscillations of the drift current 
may be induced in the conducting ring with inhomogeneous magnetic prop- 
erties under the condition tat << w"^ << (see Eq.(|6j |9|. In the case of 
the heterostructures, the inhomogeneity may be induced by the inhomogeneous 
Rashba splitting due to the space-dependent electrical field. Assuming that the 
splitting is small enough, i.e. 1\pf/pf « 1, we obtain from Eq.(|9]) that 

^ — ''^pij^ § i^^y dx]^^'^. In the case of a completely nonmagnetic and inho- 
mogeneous ring the excitation of local spin polarization also induces non-linear 
"spin-pendulum-like" oscillations, but the consideration of that goes beyond the 
framework of this paper. In the case of the open circuit, one may detect a non- 
equilibrium spin-polarization, measuring the voltage between the open ends of 
the circuit. It can be done both in the hydrodynamic and diffusion regimes. 
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